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ABSTRACT 

We report an extensive search for Lya emitters (LAEs) at z = 6.5 in the Subaru Deep Field. 
Subsequent spectroscopy with Subaru and Keck identified eight more LAEs, giving a total of 17 
spectroscopically confirmed LAEs at z = 6.5. Based on this spectroscopic sample of 17, complemented 
by a photometric sample of 58 LAEs, we have derived a more accurate Lya luminosity function of 
LAEs at z = 6.5, which reveals an apparent deficit at the bright end of ~ 0.75 mag fainter L* , 
compared with that observed at z = 5.7. The difference in the LAE luminosity functions between 
z = 5.7 and 6.5 is significant at the 3 a level, which is reduced to 2 a when cosmic variance is taken 
into account. This result may imply that the reionization of the universe has not been completed at 
z = 6.5. We found that the spatial distribution of LAEs at z = 6.5 was homogeneous over the field. 
We discuss the implications of these results for the reionization of the universe. 

Subject headings: cosmology: observation — early universe — galaxies: high-redshift — galaxies: 
formation 



1. INTRODUCTION 

The cosmic reionization was undoubtedly one of the 
major turning points in the early universe. The 
measurement of cosmic microwave background (CMB) 
temperature polarization by the Wilkinson Microwave 
Anisotropy Probe ( WMAP) implies an early reioniza- 
tion at z = 10.9±2:3 l|Page et al. 1 120061) . and the com- 
plete Gunn-Peterson trough of Sloan Digital Sky Sur- 
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vey (SDSS) QSOs suggests that c osmic reionization 
ended at z ~ 6 l|Fan et al.1 [2002). There are dis- 
putes over when and how the reionization has taken 
place, and which objects were responsible for it. Al- 
though QSOs are expected to be the main contributor 
of ionizing photons at the bright end of the luminos- 
ity function (LF) of ionizing sources, the QSO popula- 
tion alon e cannot account fo r all the required ionizing 
photons fWillott et al.ll200q) . and star- forming galax- 
ies like Lyman break galaxies (LBGs) and Lya emitters 
(LAEs) at the reionization epoch are the only alterna- 
tives that could dominate at the faint end. The cen- 
sus of observable galaxies at this ep och is sensitive to 
the ionization fraction of the univ e rse llYan fc W indhors 
200l iMalhotra fc Rhoadl l200l IBouwens et alJ I200E 



Stiavelli et alJl2005t iBunker et aljbo06j) . It is expected 
that the surrounding neutral intergalactic medium 
(IGM) attenuates the Lya photons so significantly 
that the number density decline of LAEs provides 
a useful observational constraint on the reioniza- 
tion epoch (Haiman fc SnaansllQQQt iRhoads & Malhotra 



tion e p_ _. 

l27)QjllHu"et alJl2002ft 



There have been great advances over the past three 
years in detecting distant galaxies at the edge of the 
cosmic reionizatio n era beyond z = 6 in both d r opout 
searches llDickinson et al.l 120041: IBouwens et al.l 120031 
[20051 iKneib et alJl2004HPell6 et, al.ll2004HStanwav et alJ 
2005) to find their s trong Lyman b r eaks and na rrow- 
band (NB) searches llHu et alJl200l iCubv et alJ 1 20031 
Aiiki et al.l2003tlKodaira et alJ2003tlRhoads et all2004 



Malhotra & R hoadsll200i iTaniguchi etail 120051) to de- 
tect their Lya emission lines. Complementary to these 
photometric surveys are direct spectro scopic approaches 
based on slitless spect roscopic searches llKurk e^ a l. 2004; 
Malhotr^ej^aL. 2005) and blind slit searches l)Tran et alJ 
20041: iMartin et all 12006]) . The Subaru Telescope plays 
an important role in these challenging searches for high- 
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TABLE 1 
Summary of Spectroscopic Identifications 








Observational run 


Instrument 


Ntot* N cand h Lya 


Ho 


[0 in] 


[On] 


Single Line c 


Ref. 


2002-2003 

2004 

2004 


Subaru FOCAS 
Keck II DEIMOS 
Subaru FOCAS 


20 13 10(9) d 
14 6 5 
19 3 2 



1 

1 


4° 

4 

14 




1 



3(5) d +3 c 



2 


1,2 

3 

3 



a The total number of objects for which we obtained a spectroscopic signal. 

b The total number of z6p5LAE candidates that meet our photometric selection criteria. 

c These possess neither line asymmetry with large Sw as in a LAE nor doublet features as in an [O n] emitter. 

d Although nine LAEs and five single-line objects were reported in T05, one of the single-line objects, SDF 
J132520. 4+273459, was reobserved on the 2004 Keck II DEIMOS run and was found to be a LAE. Another 
single-line object, SDF J132518. 4+272122 was concluded to be a LAE in this study based on the S w classification. 

°The spectra of four [O III] emitters and three single-line objects were obtained as for the 7VB921 strong emitter 
on the 2002 run. 
References. — (1) Kodaira et al. 2003; (2) Taniguchi et al. 2005; (3) this paper. 



z populations, in particular, in the Subaru Deep Field 
(SDF). The major goal of this project is to construct 
large samples of LBGs at z ~ 4 — 5 and of LAEs at 
z ~ 4.8, 5.7, and 6.6, and to make detailed studies 
of these very high z galaxy populations. The SDF's 
wide-field imaging increases the chance of discovering 
rare objects, such as the most distant galaxies. In ad- 
dition to the improved detectability, the wide field of 
view is less sensitive to the poten tially large cosmic scat- 
ter in the reionization history ( Barkana & Loeb 2004; 
iSomerville et al J 120041) . Following our fir st discovery of 
a couple of LAEs at z = 6.5 (z6p5LAEs; Kodair a et alJ 
120031) . iTanieuchi et all l)2005l T05) have revealed for the 
first time a statistically useful sample of nine spectro- 
scopically identified z6p5LAEs and estimated their to- 
tal a mount of star formation r ate density at this high- z 
end. iNaeao et al.1 1)2004 l2005|) also serendipitously dis- 
covered strong Lya emission at z > 6 from an i'-drop- 
selected sample in the SDF. Our SDF LAE sample was 
obtained from a general blank field without resorting to 
amplification of gravitational lensing by foreground clus- 
ters, providing reliable statistics about their number den- 
sity, LF, and cosmic star formation rate density. High- 
z surveys using gravitational lensing are complemen- 
tary to our survey, because they detect low-luminosity 
sources (lEllis et alJ l200lt IHu et all l200l ISantosI l200l 
iKneib et alJl2004|) . 

In this paper, we report the discovery of eight ad- 
ditional spectroscopically identified z6p5LAEs, which 
enables a more accurate estimation of their LF. The 
LF beyond z = 6 puts a critical constraint on the 
reionization epoch, as well as on the ionizing photon 
budget. The Lya photons are absorbed when pass- 
ing through the neutral IGM; therefore, it is naturally 
expected that the LF of LAEs should decline as it 
traces earlier times in the reionization epoch. Conse- 
quently, the observed abundance of LAEs during the 
reionization period shou ld indicate the neutral fractio n 
of IGM hydrogen x¥fi M llMiralda-Escude fc Ree3 fl998L 
iMalhotra fc RhoadJpOO'fl) and lStern et all 1)2003) found 
no significant evolution of the LF between z = 5.7 and 
6.5, implying that the neutral fraction of the universe 
is already low at z = 6.5. However, their LF estimate 
at z = 6.5 was poorly determined, since it was com- 
bined from several independent data sets with different 
selection criteria. On the other hand, there are updated 



model predict ions for the LAE's LF during the reion- 
izatio n epoch JH aiman fa Spaanslll999t [Le Dcllio u et alJ 
2005; Hai man fc Cenll2005|) . 

In addition, we evaluate the inhomogeneity of the sky 
distribution for our z6p5LAE sample. The high-z galaxy 
survey in a general field also has an advantage for de- 
termining spatial clustering. The i-dropout method gen- 
erally samples a wide redshift range at 5.7 < z < 6.2, 
which corresponds to a comoving distance as deep as 
~ 200 Mpc along the line of sight. The large-scale struc- 
ture within this large volume probed by the i-dropout 
method would be diluted in sky projection and thus can- 
not be revealed unless large spectroscopic samples are ob- 
tained. On the other hand, NB searches exploring only 
a small redshift range are more sensitive to the large- 
scale structure, although their thin slices of the universe 
are, at t he same time, likely to be affected by cosmic 
variance (Shimasa ku et all 12004'). The detection of an 
inhomogeneous distribution of galaxies beyond z = 6 
would be of great interest not only bec ause it would 
reveal the primeval large-s c ale structure (Stiavel h et alJ 
127)051 iMdhotra et a1]l2005t IQuchi et aUl2005D . hut also 
because it could provide evidence of patchy reioniza- 
tion, in which Lya flux is attenuated in one field and 
not in the other. In the reionization epoch, ionizing 
sources like L AEs would make cosmo logical H II regions 
around them l)Miralda-Escude et alJ I2000J) . The ion- 
ized spheres around adjacent LAEs will overlap, and the 
space overdensity of these ionizing sources would form 
a large enough H II region around them to allow high 
transmission of their Lya photons prior to reionization 
j|Wvithfi fc Loeh^OOalFTrrTaTietto et all200fl. The max- 
imum extent of the overlapped i onized regions is pre- 
dicted to be ~ 10 phy sical Mpc fWvithe fc Loebll2004t 
Furla rietto fc Obll2005[K which is comparable to the field 
of view of the SDF. A measurement of inhomogeneity 
in the spatial distribution of such a high-z population 
would draw a much more precise picture of the reioniza- 
tion process than has been previously possible. 

This paper is organized as follows. In § 2, we describe 
our new spectroscopically identified z6p5LAE sample. In 
§ 3, we derive the Lya LF of our z6p5LAE sample making 
comparison with LFs at z = 5.7. An estimate of inhomo- 
geneity in the sky distribution for our z6p5LAE sample 
is presented in § 4. We present the composite spectrum 
of our spectroscopically confirmed z6p5LAE sample in 
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§ 5. Some discussions of the implications for reionization 
based on our results are made in § 6, and the summary 
of the paper is given in § 7. 

Throughout the paper, we analyze in the flat ACDM 
model: fl m — 0.3, £1a = 0.7, and H = 70 /i7okms -1 
Mpc -1 . These p arameters are consi stent with recent 
CMB constraints IjSpereel et al J 12006(1 . Magnitudes are 
given in the AB system. 

2. NEW SPECTROSCOPIC CONFIRMATION 

Our z6p5LAE photometric candidate sample in the 
SDF was presented in T05, in which sample selection 
and star formation rate density were discussed. The 
sample was based on the flux excess objects in nar- 
rowband NB921 (A c = 9f96 A, FWHM=f32 A) im- 
age comp ared with the very deep broadband images of 
the SDF (Kashikawa ct al. 2004). Our comoving sur- 
vey volume was as large as 2.17 x 10 5 ft.y 3 Mpc 3 . In 
T05, we found 58 photometric candidates of z6p5LAEs 
down to A^i?921 = 26.0 (5 a) in the effective survey re- 
gion of 876 arcmin 2 ; nine of them had been confirmed 
as real by spectroscopy. In this section, we describe our 
extended spectroscopic confirmations of z6p5LAEs after 
T05. Table ^ summarizes our spectroscopic identifica- 
tions of ./Vi?921-excess objects in the SDF over the last 
three years. In summary, we have hitherto taken spec- 
troscopy for 22 objects that meet the photometric se- 
lection criteria of z6p5LAE and confirmed that 16 are 
really LAEs based on their asymmetric line profiles, one 
is an [O in] emitter, and five are faint single-line emit- 
ters. We have also included another spectroscopically 
identified z6p5LAE discovered serendipitously. 

2.1. Keck II DEIMOS Spectroscopy 

The z6p5L AE candidates we re observed with the Keck 
II DEIMOS l|Faber et all2003fl spectrograph on UT 2004 
April 23 - 24. We also allocated slits for A^B921-strong 
(z 1 — NB921 > 1) emitters, irrespective of their («' — z 1 ) 
color as a LAE criterion in order to see how our selec- 
tion criteria work. We used four multiobject spectro- 
scopic (MOS) masks with an 830 line mm -1 grating and 
a GG495 order-cut filter for each 7000-9000 s. integration 
time. The central wavelength was set to 7500 A for one 
of the four MOS masks and 8100 A for the other three 
masks. The slit width was 1".0 with 0.47 A pixel -1 , 
giving a resolving power of ~ 3600. The wavelength cov- 
erage was ~ 5000 — 10, 000 A, depending on position in 
the mask. The typical seeing size was 0".55 — 1".0 dur- 
ing the observation. Our z6p5LAEs were almost spatially 
unresolved on an A^i?921 image with 0".98 seeing size. 
Assuming that our LAEs were also spatially unresolved 
on the slits, the effective spectral r esolution may be bet - 
ter, depending on the source size l)Rhoads et al.l l2003). 
We also obtained spectra of standard stars BD +28 4211 
and Feige 110 for flux calibration. The data were re- 
duced with the spec2d pipeline 16 for DEEP2 DEIMOS 
data reduction. 

We allocated slits for 18 target z6p5LAE candidates, 
as well as 7VB921-strong emitters. Four of them were ap- 
parent [O in] emitters showing their characteristic double 

16 The data reduction pipeline was developed at University of 
California, Berkeley, with support from National Science Founda- 
tion grant AST 00-71048. 
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Fig. 1. — Spectra of eight spectroscopically confirmed z6p5LAEs. 
"F" ("D") in parentheses indicates that the object was observed 
with FOCAS (DEIMOS). The sky spectrum is overplotted on the 
bottom panel with an arbitrary flux scale. The spectrum of SDF 
J132518. 4+272122 which was identified as a LAE in this study 
already appeared in T05. 
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TABLE 2 
Spectroscopic Properties of 



6.5 LAEs 



ID a 


NAME 


z b 




/ s P cc (Lya) c 


L B 


P cc (Lya) d 


SFR s P cc (Lya) d > c 


FWHM f 


l 'll! 


OBS.s 








(10 


-18 


ergs/s/cm 2 ) 


(10 42 


'' h 70 crgs/s) 


(M /yr) 


(A) 


(km/s) 


(A) 




11 


SDF J132410.8+271928 


6.551 






5.39 




2.64 


2.40 


14.6 


477 


7.05 + 4.79 


F 


16 


SDF J132518. 4+272122 


6.547 






11.2 




5.47 


4.98 


11.0 


360 


4.75 + 2.35 


F 


18 


SDF J132520.4+273459 


6.596 






10.0 




4.97 


4.52 


8.2 


266 


11.58 + 0.54 


D 


30 


SDF J132357. 1+272448 


6.589 






10.6 




5.25 


4.78 


5.9 


192 


3.12 + 0.30 


D 


37 


SDF J132425. 4+272410 


6.552 






14.5 




7.09 


6.45 


8.9 


291 


16.34 ± 0.44 


D 


52 


SDF J132458. 5+273913 


6.541 






7.67 




3.74 


3.40 


9.7 


318 


7.63 ± 0.71 


D 


55 


SDF J132506. 4+274047 


6.546 






9.84 




4.80 


4.36 


6.3 


206 


3.77 + 0.44 


D 


59 


SDF J132432. 9+273124 


6.557 






4.50 




2.20 


2.00 


11.0 


359 


4.22 + 2.50 


F 



a The object IDs are those of T05, except ID=59, which is not listed in the photometric catalog of T05. 

b The redshift was derived from the wavelength of the flux peak in an observed spectrum assuming the rest wavelength of Lya to be 1215 A. 
These measurements could be overestimated in the case of significant damping wings by IGM. Also, the observed peak position was slightly 
shifted redward due to instrumental resolution. See Fig. |5] 

c The observed line flux corresponds to the total amount of the flux within the line profile. 

d No dust absorption correction was applied. 

c Estimated from the observed luminosities with the relation SFR(hya) = 9.1 X 10~ 43 L(Lya) Mq yr _1 as in T05. 

'Corrected for instrumental broadening. 

^Observed with FOCAS (F) or DEIMOS (D). 



lines and sometimes also H/3, one was an Ha emitter that 
has corresponding [O in] emission, four were [O n] dou- 
blets, five have apparently asymmetric single lines, and 
we did not obtain any signal from four targets. 

2.2. Subaru FOCAS Spectroscopy 

Three z6p5LAE candidates were also observed with 
the Su baru Faint Object Camera and Spectrograph (FO- 
CAS) IjKashikawa et al.ll2002|) spectrograph in the MOS 
mode on UT 2004 April 24 - 27. Although our pri- 
mary targets for this o bservation were z — 5.7 LAEs 
(Shimasak u et al] l2006) in the SDF, three slits were al- 
located for the z6p5LAE sample. We also allocated slits 
for strong NB921 emitters. The spectroscopy was made 
with a 300 line mm grating and an 058 order-cut fil- 
ter. The spectra cover 5400 — 10, 000 A, with a pixel 
resolution of 1.34 A. The 0".6 wide slit gave a spectro- 
scopic resolution of 7.1 A at 9200 A (R ~ 1300). The 
spatial resolution was 0".3 pixel -1 with 3 pixel on-chip 
binning. The integration time was 12, 000-16, 800 s. The 
sky conditions were fairly good with a seeing of 0".4-0".8. 
The data were reduced in a standard manner. We also 
obtained spectra of standard stars Hz 44 and Feige 34 
for flux calibration. 

We allocated slits for 21 targets: 14 of them were ap- 
parent [O in] emitters showing their double features and 
in some cases, H/3; one was an Ha emitter that has cor- 
responding [O in] emission; two had apparently asym- 
metric single lines; two had symmetric lines; and we did 
not obtain any signal for the two remaining targets. 

2.3. Spectroscopic Results 

We combined our spectroscopic sample with those re- 
ported in T05, in which nine LAEs and five single emit- 
ters are contained. One object, SDF J132520.4+273459, 
which was classified as a single emitter in T05, was reob- 
served with DEIMOS and was found to be a LAE based 
on the resulting better quality spectrum. We also ob- 
tained spectra for four [O ill] emitters and three single 
emitters classified as 7V_B921 strong emitters in the same 
MOS observation, although these objects do not satisfy 



the LAE criteria and were not reported in T05. The total 
spectroscopic sample for this study comprises 53 objects. 
It is difficult to identify a LAE at very high z with 
little continuum flux and a tiny signature of Lya emis- 
sion. The asymmetric line profile is the best diagnostic 
of high-z Lya emission, which results from absorption 
by neutral hydrogen; therefore, it strongly depends on 
the ionization structure in and around the high-z ob- 
jects. Althou gh some brigh t LAEs show the continuum 
breaks (Kodai ra et al.l2003fl caused by IGM attenuation, 
most are too faint to detect the break at the highest z 
epoch. We have no other spectral features but asymmet- 
ric emission profiles that can distinguish high-z LAEs 
from foreground [O n], [O in], or Ha emitters. To quan- 
tify this asymmetry accurately, we introduced an asym- 
metry statistics skewness S and weighted skewness S w . 
Here we regarded the observed spectrum, which is ba- 
sically a two-dimensional array of the flux (/j) and the 
pixel coordinate (xi), as a distribution function with an 
array size of n. The S statistic is defined as 



S 



Ia 3 ^ 



\X% Xj Ji 



(i) 



where / = J^™ /; , and W, cr are the average and disper- 
sion of Xi, respectively. The S indicator has an advan- 
tage of being independent of the line-profile modeling or 
fitting procedure. Our statistic S has a go od corre lation 
with other asymmetry indicators, a\ or a,f (Rhoads et al. 
2003J), as shown in the Appendix. 

The weighted skewness S w is the revised indicator of 
S so as to be more sensitive to an asymmetry; however, 
it does depend on the fitting procedure. We define the 
weighted skewness S w as 



5(Ai 



Al0,b)) 



(2) 



where Aio,r aud Aio,b are the wavelengths where the 
flux drops to 10% of its peak value on the red and blue 
sides of the Lya emission, respectively. 
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Fig. 2. — Thumbnail images of our eight identified LAEs. The 
object IDs are those of T05, except ID=59 which is not listed in 
the photometric catalog of T05. The B-, V-, R-, i'-, NB816-, z'-, 
and A r B921-band images are shown from left to right. Each image 
is 10" on a side. North is up and east is left. 



In this study, we classified our observed emission lines 
based on the S w indicator as shown in the Appendix. Ta- 
ble n summarizes the identifications for our 53 spectra. 
We have identified apparent foreground emitters, includ- 
ing two Ha emitters, 22 [O ill] emitters, and four [O n] 
emitters, by their multiple emission-line signatures. The 
properties of these emission-line galaxie s at z < 1.2 wi ll 
be presented in our forthcoming paper l)Lv et alJ l2006). 
The S w value of these apparent foreground emitters never 
exceeds S w = 3, which we set as the critical S w value to 
distinguish LAEs from foreground emitters. This crit- 
ical value is the same as for our z — 5.7 LAE sample 
(Shimasaku ct al. 2006). As a result, we have so far ob- 
tained 17 LAEs at z = 6.5. All of the nine LAEs identi- 
fied in T05 were classified as LAEs according to the S w 
criterion, and we obtained eight additional spectroscopic 
confirmations of LAEs in this study. 17 The spectra of 
newly identified LAEs in this study are shown in Fig- 
ure ^ an d their spectroscopic properties are summarized 
in Table |2 For all eight LAE spectra, we did not detect 
any UV continuum fluxes significant enough to measure 

17 SDF J132518. 4+272122 was classified as a single emitter in 
T05, although it shows a very red color (i' — z' > 2.21) and high 
enough S (0.173) and S w (4.75) values. We therefore classify it as 
a LAE in this study. 



Fig. 3. — Observed wavelength distribution of the Lyo line peak 
for our 17 z6p5LAE sample (solid histogram), and our foreground 
sample (dotted histogram). The transmission profile of NB921 is 
overplotted (dashed line). 



their equivalent widths reliably. Nor did we detect Ny 
A1240, the only accessible strong high-ionization metal 
line indicative of AGN activity. We discuss the compos- 
ite spectrum in § 5. Figure [5] presents postage stamp 
images of these eight LAEs in all seven bands. 

We have also eight single-line emitters in which we de- 
tected neither an asymmetric line having as large S w as a 
LAE, nor doublet features as in the [O n] emitters. Prob- 
ably most of these are unresolved [O n] doublet lines, 
based on their small S w . In fact, all of these single lines 
were observed by Subaru FOCAS, whose data lacked the 
resolving power to separate the [O n] doublet (AA = 6.64 
A) at 9160A. Our 58 photometric LAE candidates down 
to AT.B921 = 26.0 (5 a) were selected in T05 based on 
the criteria that z' - NB921 > 1, z' - NB921 > 3 a and 
i'-z' > 1.3 at i' < 27.87 (2 a) and simply z'-NB92l > 1 
at i' > 27.87 (2 a). We also adopted another criterion 
for LAE candidates with no detections (< 3 a) in deep 
B-, V-, and i?-band images. The 53 objects for which 
we obtained spectra are composed of LAE candidates 
and A^i?921-strong emitters. Twenty-two of these 53 ob- 
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jects meet the criteria for our photometric sample: 16 are 
LAEs, one is an [O III] emitter, and five are single-line 
emitters. We have another spectroscopically identified 
LAE that was not listed in the LAE candidate list by 
T05. The object is listed as No. 59 in Table and is 
found to have a very close neighbor in the i'-band image 
(see Fig. |2J), which prevented accurate aperture photom- 
etry. Sixteen LAEs meet our selection criteria out of our 
spectroscopic LAE sample of 17, indicating 16/17 = 94% 
sample completeness, provided that all eight single-line 
emitters are foreground objects. Otherwise, the sample 
completeness is (16 + 5)/(17 + 8) = 84% at most, if all of 
the single emitters are z6p5LAEs. The simple average of 
these two extreme cases, 89%, is regarded as the sample 
completeness. On the other hand, the sample contami- 
nation rate is estimated as follows. One [O ill] emitter 
satisfies our LAE criteria, 18 suggesting a 1/22 = 4.5% 
contamination rate. The contamination rate could be 
(1 + 5)/22 = 27% at most if all of these five single emit- 
ters were foreground objects. The contamination rate 
is estimated to be 16% by taking the average of these 
two cases. Therefore, our sample reliability factor, de- 
termined as the ratio of the number of true LAEs to 
the number of objects that meet our selection criteria, 
is evaluated to be (1 - 0.16)/0.89 = 0.94. One object 
found to be an apparent [O in] emitter by spectroscopy 
was removed from the photometric sample, whereas one 
LAE that happened to be found by spectroscopy but 
was not listed in the original candidate sample of T05 
was included in the photometric sample. In the follow- 
ing analysis, we used this photometric sample. There 
were six objects for which we did not obtain any signals 
in spectroscopy. Five of them have 7V_B921 > 25.5, which 
is close to the current spectroscopic limit. 

Figure |3 shows the peak wavelength distribution of 17 
confirmed LAEs, as well as 28 foreground emitters within 
the NB921 bandpass. The distribution of LAEs shows 
an apparent systematic deviation to the blue side of the 
NB transmission curve, in contrast to the symmetric dis- 
tribution of foreground emitters. This is because of the 
fact that the LAE profile, having a broad red wing and 
sharp blue cutoff, as well as a Lyman continuum break, 
makes a larger NB excess when it lies at the shorter side 
of the transmission curve. This w as also the case for LAE 
surveys at z = 5.7 with NB816 (|Shimasaku et alJl200ri 
IHu et alJl2"004|K 

3. LYa LUMINOSITY FUNCTION 

We estimated the LF of z6p5LAEs based on both our 
spectroscopic sample of 17 and our photometric sample of 
58. The LF can be simply derived from the number den- 
sity of confirmed LAEs mu ltiplied by the spectr oscopic 
confirmation fraction l|Malhotra & Rhoadsl r2004L T05); 
however, the uncertainty in this fraction inevitably de- 
pends on the magnitude because spectroscopic confirma- 
tion is more difficult at fainter magnitudes. Thus, we 
estimated the range of acceptable z6p5LAE LFs as spec- 
ified by the upper and lower limits. 

The lower limit to the z6p5LAE LF is based on our 
spectroscopic z6p5LAE sample of 17, whose Lya emis- 
sion has been securely detected, although this sample 

18 No. 29 in T05. This object is located near a bright star so 
that accurate photometry is prevented. 
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Fig. 4. — Comparison of the Lya fluxes measured by spectra, 
fspec, with those inferred from photometry, f p hot> f° r our spec- 
troscopic z6p5LAE sample. The solid line represents a one-to-one 



correspondence between f 3pec and f 



phot 



The errors were estimated 



based on the sky rms fluctuation on each spectra for f spe c, and er- 
rors on magnitudes for i p hot- One object far out of agreement is 
found to have a spectrum affected by bad columns. 



is incomplete. The upper limit was estimated from the 
combined spectroscopic and photometric samples. The 
Lya and rest-UV continuum (at 9500 A ) fluxes of our 
photometric sample were inferred from the A\B921 and 
z'-band photometry using equations (7) and (6) of T05, 
respectively. Note that we derived the apparent Lya 
luminosity uncorrected for either the dust extinction or 
the self-absorption evident on the blue-side cutoff of the 
emission line. The comparison of Lya fluxes measured 
spectroscopically with those inferred from photometry 
for our spectroscopic z6p5LAE sample is shown in Fig- 
urc0] The correspondence is good except in a few cases. 
One object far out of agreement is found to have a 
spectrum affected by bad columns. We therefore used 
the photometric inferred Lya fluxes for the remaining 
41 (= 58 — 17) objects in the photometric sample. 

We have to correct for detection completeness, which 
could crucially affect the result when calculating the LF 
based on a deep photometric catalog. The detection 
completeness as a function of apparent NB921 mag- 
nitud e was estimated in almost the same way as in 
Kashikaw a et alj <|2004j) , that is, by counting detected ar- 
tificial objects distributed on the real NB921 image. We 
assumed Gaussian profiles of FWHM= 1".13 for these ar- 
tificial objects, which is the nominal size of our z6p5LAE 
sample objects (T05). The detection completeness was 
thus evaluated as > 75% for NB921 < 25.0 and 45% at 
the limiting magnitude of NB921 — 26.0. In the upper 
limit estimate, we corrected for this detection complete- 
ness by number weighting according to the NB921 mag- 
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Fig. 5. — Cumulative Lyo LF of our z6p5LAE sample. The 
open circles denote the raw counts of our spectroscopic sample + 
additional photometric sample, and the filled circles are corrected 
for detection completeness (upper limit). The triangles denote the 
raw counts of the pure spectroscopic sample (lower limit). Error 
bars for the filled circles are Poissonian. The squares and long- 
dashed li ne indicate the LF of L AEs at z = 5.7 evaluated from 
the SDF IShimasaku et al.l f2006V The short-dashed and dotted 
lines show the Schechter LFs, in which the Lya luminosities are 
reduced by a factor of 0.6 (L* X 0.6) and 0.4 from the z = 5.7 LF, 
respectively. [See the electronic edition of the Journal for a color 
version of this figure^ 



nitude. Note that the upper limit can be regarded as our 
current best estimate for the z6p5LAE LF, because it is 
properly corrected for the detection incompleteness. 

Figure[S]shows the cumulative Lya LF of our z6p5LAE 
sample. The open circles denote the raw counts of the 
spectroscopic + additional photometric sample, and the 
filled circles are those with corrected detection complete- 
ness. The triangles denote the raw counts of the pure 
spectroscopic sample. Therefore, the triangles and cir- 
cles are the lower and upper limits, respectively, of our 
estimates of the z6p5LAE LF. The error bars on the 
filled circles just denote the Poisson errors, although 
there may be other plausible error sources, such as an 
ambiguity in inferring the Lya luminosity from photo- 
metric data. Taking into account the corrections with 
respect to the sample reliability factor (94%) evaluated 
in the previous section, the LF has a margin to go up- 
ward by a factor of 1.06, although this uncertainty is 
smaller than the Poisson errors. We fitted a Schechter 
function, <j)(L)dL = (j)*{L/L*) a exp(-L/L*)dL/L, to our 
z6p5LAE LF. The % 2 was minimized with a single grid 
search in the two parameters, L* and <f>* , for fixed slopes 
of a = —2.0, —1.5, and —1.0. The derived best-fit pa- 



rameters are listed in Table [3J 

The squares and dashed line in FigurcEl show the LF of 
LAEs at z = 5.7 evaluated in the SDF IjShimasaku et all 
2006). This is the most reliable estimate so far be- 
cause of the deeper photometric sample and larger num- 
ber of spectroscopic confirmations than that of lHu et alJ 
(2004). In addition, it is worth noting that this LAE 
sample at z = 5.7 was selected from the same field and 
the same photometry as this study, having similar sur- 
vey volume and similar selection criteria. For example, 
we have carefully determined the NB-excess criteria so as 
to have almost the same equivalent width (EW> 20 A in 
rest frame) thresholds for both of these LAE samples. 

According to previous studies, the LF o f LAEs be- 
tween z = 3 and 6 is almost unchanged (Tr an et alJ 
120041 Ivan Breukelen et~aT]|2005l) . The LAE population 
appears to have a quite similar number density at all 
epochs. A comparison of the LAE LF from previous 
studies at z < 5 . 7 is s hown in a companion paper by 

which confirmed that the LF of 
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LAEs at z — 5.7 is almost identical to those in lower 
z. However, our z6p5LAE LF shows an apparent deficit 
at the bright end in both upper and lower limit esti- 
mates. The L* difference between z = 6.5 (the upper 
limit) and z — 5.7 is a factor of ~ 2, corresponding to 
~ 0.75 mag, assuming a fixed a = —1.5. We have al- 
ready spectroscopically identified almost all of our bright 
LAE candidates (four confirmed out of six candidates at 
A^i?921 < 25.0), and so this resulting deficit at the bright 
end would not change significantly, even if we obtain 
more follow-up spectroscopy in the future. Moreover, un- 
certainties in the detection completeness correction and 
the selection effects according to the equivalent width of 
the Lya line are expected to be small at the bright end. 
As seen in Figure 0J there are possible but unsystematic 
errors in inferring the Lya luminosity from photometric 
data, which could affect the resulting LF to some degree. 
When using only the photometric inferred Lya luminosi- 
ties for our entire 58 object sample, the best-fit Schechter 
parameters only change by at most Alog(i*) = 0.04 and 
Alog(</>*) = 0.1 for -2.0 < a < -1.0, which is negligible. 
Moreover, we carried out a Monte Carlo simulation to 
investigate any possible distortion that the discrepancy 
between spectroscopically measured and photometrically 
inferred Lya luminosities could cause in the resulting LF. 
We re-calculated the LF many times after assigning a 
random error perturbed within the same scatter as in 
Figure 01 to each Lya luminosity. With this simulation, 
the best-fit Schechter parameters vary with rms fluctu- 
ations of cr(log(£*)) = 0.019 and a (log (<£*)) = 0.032 for 
fixed a = —1.5, suggesting again that the errors in infer- 
ring the Lya luminosity are expected to be small. 

At lower z (z = 3.0 — 5.7), several LAEs with large 
L(Lya)> 2 x 10 43 h^ ergs s _1 have been actually iden- 
tified bv_srje£tros£Ojjy^with much higher number den- 
sity (jS himasa ku et aljl2006t IHu et alJl20Q4 iMaier et alJ 
l2003HKudritzki et al.ll2000HCowie h Hulll998fl . whereas 
our z6p5LAE sample includes no spectroscopically con- 
firmed objects of such large Lya luminosity. Three other 
LAEs at z ~ 6.5 have been spec troscopically i denti- 
fied so far in independent surveys JRhoads et al. 2004; 
iKurket al.ll2004t iStern et al.ll2005l) without taking ad- 
vantage of the amplification by a foreground gravita- 
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TABLE 3 
Best-fit Schechter Parameters for LAE LF at 
5.7 of the SDF 



: 6.5 AND 



Sample Limit 



z = 6.5 Upper 



Lower 



a L* cj>* 

(fix) log(//i" 2 ergs s" 1 ) log(/ft|„ Mpc" 3 

or. Ar, „j+0.14 n 1 „+0.30 



-2.0 
-1.5 
-1.0 
-2.0 
-1.5 
-1.0 



42.74 
42.60 

42.48 



+0.14 
0.14 

+0.12 
0.10 

+0.10 
0.08 



43.30±g;2S 

42.95±g;™ 

42.75ig;j 



3.14 

-2 

-2.74 

-4.80 

-4.17 

-3. 



5.7 



-2.0 
-1.5 
-1.0 



43-30t£^ 
43.04±H| 

,+0.10 
*-0.10 



+0.30 

0.34 
+0.24 
-0.26 
+0.18 

0.22 
+ 1.02 

0.20 
+0.70 

0.83 
+0.51 

0.57 



12.1 



-3.96 
-3.44 
-3.14 



+0.28 
0.30 

+0.20 
0.16 

+0.14 
0.12 



tional lens. All of these LAEs also have similar Lyo 
luminosities of L{Lya)= 1.04 — 1.1 x 10 43 h^ ergs s -1 , 
consistent with our brightest LAEs. 

On the other hand, at fainter luminosities of L(Lya)< 
5 x 10 42 hjQ ergs s _1 , our upper limit estimate of the 
z6p5LAE LF almost reaches the same amplitudes as that 
at z = 5.7, although the amplitude difference between 
our upper and lower limits is too large to constrain its 
faint end. Our faint spectroscopic sample is still too small 
to establish whether there is a significant difference be- 
tween the faint end of the LF at z = 6.5 and 5.7. In 
this study, we cannot determine the faint-end slope of 
the LAE LF at z = 6.5, and consequently it is difficult 
to constrain the true contribution of the LAE population 
t o the entire photon bu d get re quired for full reionization. 

iMalhotra fc Rhoadsl 12004]) and iStern et alJ pOOl 
found no significant evolution of the LAE LF from 
z = 6.5 to 5.7. However, their estimates were based 
on small samples combined from various independent 
data sets with different selection criteria. Our z6p5LAE 
sample was selected with the same criteria from a large 
homogeneous sample in a general field without resort- 
ing to amplification of gravitational lenses. The uncer- 
tainty of our detection completeness estimate is small, 
at least at the bright end. Moreover, our survey comov- 
ing volume is c ompara bly as large as o thers at z = 5.7 
(Shimasa ku et alJ l2006: Hu et alfeoO^ . As seen in Fig- 
ure |31 our effective survey depth is smaller than that 
estimated from the FWHM of the NB filter; however, 
it is more or less the same as for the z — 5.7 LAE 
samp le, and its correc tion would not affect the result 
(Shimasaku ct al. 2006). Therefore, our LF estimate of 
z6p5LAEs is highly reliable, although we cannot com- 
pletely rule out the possibility that the deficiency at the 
bright end is caused by cosmic variance. 

To illustrate the significance of the LF difference be- 
tween z = 6.5 and 5.7, we plot the error contours for 
our Schechter-parameter fits in Figure [H| In this case, 
we compare only the upper limit LF estimate, which 
is the current best estimate for our z6p5LAE sample 
and is most appropriate for establishing the significance 
of the LF difference compared to z = 5.7. The confi- 
dence levels of the fitting were computed based on Pois- 
sonian error statistics. The best-fit parameters of the 
z = 5.7 LF are slightly different from those presented in 



43.5 



43 



42.5 





43.5 


* 


43 


Q0 


42.5 



43.5 



43 



42.5 




= 6.5 LAET 



"a = -2.0 

I I I I I I I I 



■+— H 





7 a = -1.0 



Fig. 6. — Error ellipses of the best-fit Schechter parameters 
4>* and L* for each fixed o-value. The lower set of ellipses in 
each panel are for the z6p5LAE upper limit sample, while the up- 

§ er set of ellipses in e ach panel are for the z5p7LAE sample of 
himas aku ct al. (2006). The inner and outer solid ellipses are the 
1 a and 3 a confidence levels, respectively. The dotted ellipses 
show the 2 a confidence levels. The error bars in the lower panel 
(left error bar, z5p7LAE sample; right error bar, z6p5LAE sam- 

t le) correspond to the uncertainty in <j>* due to cosmic variance 
somervillc ct al. 2004). [See the electronic edition of the Journal 
for a color version of this figure.] 



Shimasak u et al.l (|2006), in which Schechter parameters 
were determined so as to be consistent with the NB816 
number count. In contrast, here we determined these 
parameters simply by fitting a Schechter function to the 
data points. Figure reveals that the (L*, <fi*) error el- 
lipses for z = 5.7 and 6.5 do not overlap each other for 
any a; that is, the difference in LF between z = 5.7 and 
6.5 is significant at more than the 3 a level. The diffcr- 
en ce in L* is more s i gnific ant than that in the <fi* . Based 
on lSomerville et all I 2004J) . we evaluated the cosmic vari- 
ance of our z6p5LAE sample. We assumed a one-to-one 
correspondence between LAEs and dark haloes, and used 
their predictions at Z = 6. With the comoving survey 
volume of 2.17 x 10 5 ft.f 3 Mpc 3 and the number density 
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of 2.67 x 10" 



'70 



Mpc~ J (7.83 x 10"°) for the upper 



(lower) limit estimate, we obtained a cosmic variance of 
~ 32% (~ 37%). We also estimated a variance of - 20% 
for the z = 5.7 LAE sample. The possible field-to-field 
varia nce in the LF at z = 5 .7 can be seen in Figure 
11 of Shimasaku et al. (2006). As shown by the error 
bars in Figure El the 3 a error circles for the two epochs 
overlap each other when taking into account the cosmic 
variance; however, our upper limit estimate still differs 
from the z = 5.7 result at the 2 a level. 

4. REST-UV CONTINUUM LUMINOSITY FUNCTION 

In the previous section, the flux of the rest-UV con- 
tinuum (at 9500 A) was simultaneously derived from the 
iVi?921 and z'-band photometry. We derived the rest- 
UV (1255 A at z — 6.57) continuum LF based on our 
photometric sample of 58 objects. The correction for de- 
tection incompleteness in 7V.B921 was taken into account 
when calculating the LF, although the correction should 
actually be based on the completeness measured in z' 
band, which corresponds to the rest-UV flux. However, 
as our Lya-selected sample is basically constructed from 
an 7V_B921 magnitude-limited sample, it is inevitable 
that the derived rest-UV continuum LF may be affected, 
especially at the faint end of the LF, by the difference in 
completeness of the iVi?921 and z' bands. This is also 
the case for ot her LAE studies (Shimas aku et al]l2006t 
IHii fit. a,Ul200l. 

Figure shows the rest-UV continuum LF of our 
z6p5LAE sample compared with other studies at similar 
redshifts. No correction has been applied for dust. The 
vertical lines indicate the corresponding limiting mag- 
nitudes in z' band. Our LF measurements at magni- 
tudes fainter than Mjjy — —20.24 (3 a) may be uncer- 
tain because the corresponding z'-band magnitudes are 
no longer reliable. We overplot in Figure other rest- 
UV conti nuum LF estimates of t he SDF LAE sample at 
z = 5.7 (Shimasa ku et al J 120061) . the i-dropout objects 
at z ~ 6 of iBouwens et alJ 1 2005 ). and the LAE sam- 
ple at z = 5.7 of IHu et all 1 2004 ). We neglect here a 
slight difference in corresponding rest-frame wavelengths 
(~ 1350 A at z = 5.7 and - 1255 A at z = 6.5), as- 
suming a flat far-UV spectral energy distribution. Our 
measurements agrees with these two studies very well at 
Muv < —20.5. The agreement of the rest-UV contin- 
uum LF at the bright end for z — 6.5 and 5.7 is in clear 
contrast to the difference seen in the Lya LF. It should 
be noted that the rest-UV continuum luminosity is not 
attenuated by the neutral IGM and is less affected by 
dust extinction than the Lya luminosity. As far as the 
rest-UV continuum LF is concerned, cosmic variance is 
not severe for these samples. The flatter faint-end slope 
of the LFs of LAE samples at both z = 6.5 and 5.7 com- 
pared to that of i-dropouts at z ~ 6 could be caused by 
the detection incompleteness of the LAE sample. 

5. CLUSTERING PROPERTIES 

Our sample of 58 z6p5LAE candidates have been ex- 
tracted from a very wide field of view (34'x27'). We tried 
to detect a clustering signal in the z6p5LAE sample using 
several methods. We derived the angul ar two-point cor- 
relation function (ACF) w{9) using the Landv & Szalay 
( 1993) estimator. One hundred thousand random points 
were created with exactly the same boundary conditions 
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Fig. 7. — Rest-UV (1255 A at z = 6.57) continuum LF of our 
z6p5LAE sample (circles and solid line). As a comparison, the 
squares (and short- dashed line) are the rest-UV LF of the LAE 
sample at z = 5.7 evaluated in the SDF (Shimasaku et al,l l2006fl . 
the trian gles (and long- dashed line) are the «-dropout objects at 
z ~ 6 of Bouwc ns et al.l 120051) . and the crosses (and dot -dashed 
line) are the LAE sample at z = 5.7 of IHu et aT] (2004). The 
vertical lines indicate the limiting magnitudes in z' band at Mjjy = 
-19.05, -19.80 and -20.24 for 1 a, 2 a, and 3 <r, respectively. Error 
bars are Poissonian. [See the electronic edition of the Journal for 
a color version of this figure.] 



as the SDF galaxy catalog, avoiding the mask regions in 
which saturated stars dominate. The top panel of Fig- 
ure [S] shows the ACF for z6p5LAE. Circles denote the 
ACF for the 58 objects of z6p5LAE, whereas squares 
denote the 53 objects of the z6p5LAE sample, exclud- 
ing five single-line emitters. We did not correct for 
the integral constraint, which is negligible in the SDF 
(Kashika wa et al.ll20 06). We estimated only Poissonian 
errors on the ACF as a w {9) = {(1 + w(9))/DD(6)} - 5 , 
where DD(9) is the number of random-random pairs hav- 
ing angular separation 6. The result shows that the am- 
plitude is almost zero for all scales, indicating that the 
sample has an almost homogeneous distribution. How- 
ever, our z6p5LAE sample is so small that the derived 
ACF has a large ambiguity. Therefore, we also applied 
two other methods that are more robust for small num- 
ber statistics to quantify the clustering strength. 

First, we applied the two-dimensional Kolmogorov- 
S mirnov test to o ur sample. This test was generalized 
bv lPeacockl l)1983|) to give the integral probability distri- 
bution in four quadrants around a certain point. To see 
the difference from a homogeneous distribution, we gen- 
erated random points as for our ACF estimate. We found 
that our z6p5LAE sample was equivalent to a homoge- 
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Fig. 8. — Top: ACF of our z6p5LAE sample. The circles denote 
the ACF for the entire z6p5LAE sample of 58 objects, whereas 
squares denote that for the 53 z6p5LAE sample in which five single- 
line emitters were removed. Error bars show the 1 <r Poissonian 
errors. Bottom: VPF of our z6p5LAE sample (circles). The dashed 
line is the VPF of a random sample. 



neous distribution at the 83.3% (83.8% after removing 
the five single-line emitters) confidence level. 

Second, we estimated the void probability function 
(VPF). The VPF is defined as the probability of hav- 
ing no galaxies in a randomly placed sphere of radius 
R, or in a circle of angular radius 9 in the case of a 
two-dimensional sky distribution. The VPF is known to 
be rel ated to the h ierarchy of n-point correlation func- 
tions j^y hit^ejll979'l . We adopted the same technique as 
iCroton et al.l l)2004|) to correct for the irregular geometry 
of the survey region. The result is shown as a solid line 
in the bottom panel of Figure |H1 compared with that of a 
random sample indicated by the dashed line. No excess 
in the VPF was found for our z6p5LAE sample relative 
to the random sample. 

These three estimates indicate that the spatial distri- 
bution of the z6p5LAE sample is homogeneous. The 
results were identical for all three estimates, even if 
we divided our sample into brighter/fainter subsamples. 
The homogeneous distribution of z6p5LAEs is in con- 
trast to that at lower z, where LAE populations are 
often found to trace the large- s cale structures even in 
blank fields (Igteidel et alJ 1200(1 jShimasaku etTaf] 120031: 
iPalunas et al.ll2004HOuchi et alJl2005|) . 

As in the previous section, we estimated the cosmic 
variance as ~ 32% at most. We should note that NB 
searches exploring only a small redshift coverage are also 
sensitive to the large-scale structure (Shimasaku et al. 



|2004[) . We cannot rule out the possibility that we happen 
to see a very homogeneous region of the z = 6.5 universe. 
A survey of a much larger volume is required for further 
discussion. 

6. LYa PROFILE OF THE COMPOSITE SPECTRUM 

Although each individual spectrum of our 17 
z6p5LAEs has too low a signal-to-noise ratio (S/N) to 
infer either the internal dyna mics of the LA E itself or 
the IGM characteristics (e.g.. lHaiman 1 12002'). the com- 
posite spectrum could be useful to see the general spec- 
troscopic properties of the z6p5LAE population. We 
have 17 z6p5LAE spectra of different spectroscopic reso- 
lutions. First, we removed five spectra 19 that have been 
taken with the poorest instrumental resolution. Each 
spectrum was then smoothed with an adequate Gaus- 
sian kernel chosen to produce a common instrumental 
resolution of FWHM= 6.41 A. The instrumental reso- 
lution for each spectrum was practically measured from 
the FWHM of sky lines near the Lya emission. The red- 
shift was measured based on the line peak wavelength, 
then shifted to the mean redshift (z) = 6.556, and re- 
binned to a common pixel scale. In the process, we 
neglected possible systematic offsets of Lya lines from 
the rest frame establis hed by other lines, w hich are often 
found in LBG spectra IjShaolev et al.l20f)3t) . The spectra 
were then combined by taking the average after scaling 
according to their peak line flux. 

The top panel of Figure OH shows the final composite 
spectrum from our 12 z6p5LAE sample. The composite 
spectrum reveals an apparently asymmetric profile with 
an extended red wing, which is shown clearly by a com- 
parison with a Gaussian profile (dotted line) correspond- 
ing to the final instrumental resolution. The skewness 
and weighted skewness of the composite spectrum were 
S = 0.542±0.007, and S w = 11.466±0.156, respectively. 
The blue side of the line profile is almost adequately ex- 
plained by th e inst rumental resolu tion blur, as was als o 
concluded bv lHu et alJ l|2004fl and iWestra et alJ (2005). 
Assuming that an intrinsic Lya profile is a simple Gaus- 
sian at almost the same peak position as observed, and 
completely truncated at the blue side of the line, the 
resulting profile convolved with the instrumental resolu- 
tion did not coincide with the observed profile, instead 
producing a large red-wing anomaly. This disagreement 
is inconsistent with the results for the z = 5.7 case by 
IHu et al] pOOl . 

To explain the observed profile, we considered two 
plausible models. The first model is the "galactic wind 
model" which was motivat ed by the similar analogy with 
Dawson et all l)2002|) and iMas-Hesse et all l|2003|) . If 
galactic winds are present, the far side of the expanding 
shell back-scatters redshifting Lya photons that would 
make another broadly extended Gaussian component in 
their line profile. We simply assumed that the Lya pho- 
tons from the blue side of the object redshift are com- 
pletely absorbed by neutral hydrogen at the near side in 
lines of sight. This model is composed of two Gaussian 
profiles: one is a high-amplitude narrow Gaussian that 
originates from recombination Lya photons in the cen- 
tral H II region, and the other is a low-amplitude broad 
Gaussian from Lya photons back-scattered by a galactic 

19 The removed sample were object ID=4, 5, 6, 7, and 8 in T05 
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TABLE 4 
Model Parameters for Composite Lyo Profile of LF of z6p5LAE 



Model 


Component 


(A) 


Amplitude 
(ergs s — 1 cm -2 A -1 ) 


(A) (kms" 1 ) 


Offset 
(A) (km s" 1 ) 


Galactic wind model 


Central Hn 
Galactic wind 


9183.3 
9191.6 


2.11E-17 
3.31E-18 


3.32 108.5 
7.93 259.2 


-2.34 -76.5 
+5.97 +195.1 


Reionization model 




9179.9 


9.00E-17 


6.50 212.4 


-5.70 -186.3 
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Fig. 9. — The composite spectrum of our 12 z6p5LAE sample. 
The spectrum in the top panel shows the final composite spectrum. 
The dotted line denotes a Gaussian profile corresponding to the in- 
strumental resolution. The middle panel shows the best-fit galactic 
wind model. The resulting profile is shown as a smooth solid line 
that is the combination of narrow and broad Gaussians convolved 
with the instrumental resolution. The lower panel shows the reion- 
ization model. The resulting profile is shown by a smooth solid 
line that is the combination of the intrinsic Lya profile and damp- 
ing wing attenuation convolved with the instrumental resolution. 
See Table |41 for the best-fit parameters for each model. [See the 
electronic edition of the Journal for a color version of this figure.] 



wind. The middle panel of Figure ED shows the best-fit 
galactic wind model. The resulting profile is shown by 
the solid line, which is the combination of narrow and 
broad Gaussians convolved with the instrumental reso- 
lution. It perfectly explains the observed profile. The 
best-fit parameters are listed in Table 0] The picture of 
large-scale outflowing of gas with velocities ~ 200km s _1 
is in good agreeme nt with those of nearby H II galaxies 
l|Kunt,h et al.lll998fl but is somewhat smaller than those 
of z ~ 3 LBGs llPettiniet alJ l2002fl an d z ~ 5 LAEs 
Ipawson et al.ll2002HWestra et alJl2005|) . 
The second model is the "reionization model" , in which 



the intrinsic Lya line has a larger amplitude than that 
observed and its peak wavelength is much shorter than 
the observed peak position. Such a picture is generally 

S iredicted for the Lya profile in the reionization epoch 
Haiman 1120021 ISantosll2004fl . The Lya photons would 
be absorbed by both the red damping wing of the Gunn- 
Peterson trough from outside the cosmological H II region 
and the residual neutral hydrogen inside the H n region. 
The latter's true contribution is not yet known; thus, we 
simply assumed that the inside of the H n region is suf- 
ficiently optically thick to completely attenuate the blue 
side of the line. The damping wing scattering would be 
effective at wavelengths larger than that corresponding 
to the blue edge of the H n region; thus, the red side of 
the line would be attenuate d by the damping;. We used 
the damping optical depth of lLoeb et alJ (2005). We fur- 
ther assume d the radius of the H II reg ion to be 0.45 
proper Mpc (|Haiman 1 2002: Santos 200j|. 

The bottom panel of Figure El shows the cosmological 
H II region model fitting to the data, and the best-fit 
parameters are listed in Table 0] This model also ex- 
plains the observed extended red wing fairly well. The 
predicted intrinsic Lya luminosity from this model is 4.6 
times that observed, which is roughly consistent with 
the case of IHu et alJ |2002) and its model prediction 
bv lHaiman I <|2002) . However, the transmission factor 
of ~ 20% of the total line flux is smaller than the factor 
of ~ 50% dimming suggested by the LF. Here, we neglect 
the possible luminosity dependence of Lya attenuation, 
which was suggested by the LF shown in § 3. The dis- 
crepancy in Lya attenuation between that implied by the 
line-profile model fitting and that from the LF difference 
can be reduced by assuming a larger radius for the H n re- 
gions. For example, we can obtain a ~ 40% transmission 
factor if the radius of the H n region is as large as ~ 0.90 
proper Mpc. Although we assume an optically thick core 
for the H II region, any escape of Lya flux at the blue 
side would broaden the line profile, making the fit worse. 
We conclude that a larger contribution of the residual 
neutral hydrogen inside the H n region compared to the 
damping wing is required, although our data still lack 
the spectral resolution to make quantitative predictions 
of the density profile inside the H n region. 

A better fit is obtained for the galactic wind model; 
however, our composite spectrum still has too low S/N 
and spectral resolution to determine which model is more 
plausible. It should also be noted that there is likely to 
be a scatter in FWHM among our LAE sample, and so 
it is unclear whether all of the LAEs have prominent red 
wings that appeared in the composite spectrum. While 
in the model, the radiative transfer process of Lya pho- 
tons through the emitting galaxy and the IGM is too 
complicated to justify every parameter based only on the 
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line-profile fitting. 

7. IMPLICATIONS FOR REIONIZATION 

In this study, we found that the z6p5LAE LF has a 
clear deficit at its bright end compared with that at z = 
5.7. The simplest interpretation is that the LAE popu- 
lation undergoes Lya luminosity evolution from z = 6.5 
to 5.7. Strictly speaking, the LAE population has some 
evolution in its EW from z = 6.5 to 5.7, given that we 
observe no evolution in the rest-UV continuum LF. The 
LAEs must be a very young population, having ongo- 
ing starbursts in so short period that even 100 Myr is 
a long time, over which their Lya luminosity can easily 
drop. Nevertheless, the number density of LAEs does 
not change from z ~ 3 to ~ 5.7. Thus, it is more natural 
to assume that this lack of LF evolution should extend 
up to z — 6.5, as opposed to the LAE population having 
strong evolution between z — 5.7 and 6.5. The num- 
ber density decline from z = 5.7 to 6.5 could imply a 
substantial transition in the cosmic ionization state be- 
tween these epochs. In this section, we offer a possible 
interpretation of these observational results in the con- 
text of the reionization of the universe, assuming that 
the nature of LAEs themselves has no drastic evolution 
between z = 5.7 and 6.5. 

Assuming a fully ionized IGM at z = 5.7, the compar- 
ison of the LFs at z = 5.7 and 6.5 puts constraints on 
the neutral fraction of IGM hydrogen x^ M at z = 6.5. 
The short-dashed and dotted lines in Figure [5] show the 
Schechter LF, in which the Lya luminosities are reduced 
by a factor of 0.6 (L* x 0.6) and 0.4 from the LF at 
z = 5.7, which is still consistent with our upper and 
lower limit LF estimates, respect i vely. According to the 
IGM dynamical model of lSantoa 1)20041 ) . a Lya luminos- 
ity drop by AL* - 0.75 from the fully ionized IGM 
corresponds to x^ M = 0.45. However, this predicted 
v alue is s t rongly model dependent, and even the model 
of lSantoa 1)20041) covers a wide range of acceptable mod- 
els with different predictions of Lya attenuation. The 
predicted value of x^ M can be much smaller than 0.45 
in some of these models. Therefore, our LF estimate 
could allow a neutral fraction of the IGM at z = 6.5 of 
< Zm M ;$ 0.45. This upper limit of x]^ M at z ~ 6.5 is 
consis tent with the recent resu lts of lMalhotra fc Rhoada 
(2006 *) andlTotani et all l)2006fl . 

Hai man &: Ce n (2005) evaluated the LF evolution of 
LAEs, taking into account the luminosity dependence 
of the Lya flux attenuation. In the epoch of reioniza- 
tion, ionizing sources like L AEs would make surroundin g 
cosmological H n regions (Miralda- Escude et alJ 120001. 
Based on the CDM model, the galaxies embedded in mas- 
sive dark halos would collapse first, so that the ionizing 
sources in this era would be preferentially located in the 
high-density regions. The H n regions of the bright LAEs 
clustered in the overdense regions would overlap effec- 
tively and create a larger H n region with a high ioniza- 
tion fraction, which would significantly reduce the Lya 
flux attenuation. As a result, it is predicted that bright 
LAEs should be readily observed, whereas faint LAEs 
are more severely attenuated. However, this luminosity 
dependence would just shift the Schechter LF downward 
by a certain factor according to xj^f M because the LF is 
steeper at the bright end. Our observed spectroscopic LF 
(lower estimate) is nearly consistent with such a trend, 



but our upper limit estimate of the z6p5LAE LF, which 
has a steep decline only at its bright end, is not. Our 
null result of finding any signals of bright LAE clustering, 
is also inconsistent with this model. As Haiman & Cen 
(2005) suggested, the predicted LF profile strongly de- 
pends on the model assumptions, such as the constant 
escape fraction with respect to Lya luminosity. 

Assuming a clumpy IGM and discrete ionizing sources, 
there are two conflicting model predictions about the 
spatial inhomogeneity of reionization propagation. In a 
first phase of reionization, most of the sources are formed 
in high-density regions and ionize the dense gas around 
them. If the local neutral IGM around these sources is 
dense enough to allow higher recombination rates than 
ionization rates, the reionization will not complete in 
these high-de nsity regions. Consequently, a s in the model 
prediction bv iMiralda-Escude et al.l (2000), H II regions 
would expand preferentially towards low gas density re- 
gions, and overdense regions would gradually be ionized 
after the epoch of overlap of cosmological H II regions 
( "outside-in" model). This picture is supp orted by hy- 
drody namic IjGnedinl |2000) and A-body l|Ciardi et alJ 
2003) simulations. On the other hand, if the local IGM 
density is low enough to help an efficient escape of Lya 
photons, it is expected that the dense IGM region, in 
which a large H II region can be formed by overlap of ion- 
ized bubbles around high-luminosity sources, would be 
ionized first. It would then proceed to a void where it is 
dominated by only low-luminosity sourc es ( "inside-out" 
model). S uch a picture is p redicted by Sokasian et alJ 
(2003J) and lFurlanetto et alJ 1)2006)) . 

The key diverging point between these two 
models is the IGM den sity distribution, which 
IMiralda-Escude et alJ 1)2000)) assume as an extrapo- 
lation from that at z = 3; howe ver, the adequacy of that 
assumption is still uncertain (Furlanetto et al. 2006). 
These two contradictory models obviously have different 
predictions about the spatial distribution of observable 
galaxies at that epoch: galaxies residing in the under- 
dense regions are easily observed in the outside-in model, 
whereas galaxies residing in the overdense regions are 
easily observed in the inside-out model. In turn, a more 
inhomogeneous galaxy distribution is predicted in the 
inside-out model. Assuming that we are seeing the final 
stage of reionization through our observed z6p5LAE 
sample as deduced from the previous estimation of 
< x^f M < 0.45, the homogeneous distribution found in 
this work agrees better with the outside-in model. Our 
assumed picture, in which the LAE overdense regions 
were still obscured by surrounding thick neutral IGM 
clouds at z — 6.5, is also consistent with our LF deficit 
in bright LAEs. Such an overdense region could be 
kept neutral until later epochs, trapping bright ionizing 
sources like luminous QSOs, in which the appearance 
of Gunn-Pet erson troughs are detected at lower z : 
z ~ 5.2-5.8 l|Becker et al.ll2001llDio7govski et all2 001). 
Inhomogeneous reionization is also suggested by the 
significant variation in the I GM transmission among 
different QSO lines of sight IjDiorgovski et al.l 12006) . 
In summary, implications for the reionization process 
suggested by this study are a high clumping factor of 
the IGM, and inhomogeneous reionization. Although 
we do not draw any further quantitative conclusions, we 
may be looking at only low-luminosity LAEs residing in 
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Fig. 10. — Correlation between the Lya luminosity and the local 
surface number density for our 58 z6p5LAE sample. The circles 
are from the spectroscopic sample, whereas the squares are from 
the photometric sample. 



Fig. 11. — Correlations between the Lya luminosity and the 
FWHM of the Lya emission line for our sample of 17 spectroscop- 
ically confirmed z6p5LAEs. The dotted line is the fiducial linear 

fit. 



low-density IGM regions at the end of the reionization 
epoch. 

We plot in FigureEHthe relation between the local sur- 
face number density of LAEs and their Lya luminosities, 
L(Lya). The local surface number density is measured 
by the number of LAEs in a circle of 8/ifg 1 Mpc comoving 
radius around each sample object. There is no apparent 
correlation between the local density and L(Lya), again 
suggesting spatial homogeneity. It should be noted that 
the local number density in Figure E3 just accounts for 
our LAE sample, and at this point, there is no way of 
inferring the presence of other ionizing sources like i'- 
dropout galaxies without strong Lya emission. Nor did 
we find any clear relation between the local densities and 
the FWHM of Lya lines. 

FigurelTTIshows the FWHM of Lya emission as a func- 
tion of L(Lya) for our spectroscopic z6p5LAE sample. 
There is a weak anti-co r relati on between FWHM and 
L(Lya). lHaiman &; Cenl l)2005|) predicted a similar cor- 
relation, which they proposed as an independent diag- 
nostic of x}§ M , aside from the LF. The blue side of Lya 
emission is dominantly attenuated by the residual H I in- 
side the cosmic H n regions. A systematic high velocity, 
as in galactic winds, would reduce more effectively the 
amount of residual H I for smaller H II regions surround- 
ing low-luminosity sources; thus, the anti-correlation is 
expected. This effect would depend more sensitively on 
the line width at the blue side, which is difficult to mea- 
sure precisely. The correlation is expected to be steeper 
for higher a;^f M ; however, the observed correlation in 



Figure ITT1 shows too much scatter to determine £^f M . 
The observed anti-correlation could be a sign of a high 
x^f M , indicating that cosmic reionization has not been 
completed at z = 6.5. A possible variation in systematic 
internal velocity among LAEs would also dilute the rela- 
tion. The intrinsic relation between FWHM and L(Lya) 
has not been clearly established. Matsuda (2005) found 
a clear positive correlation between the velocity disper- 
sion and L(Lya) for their large Lya blob sample at z ~ 3 
in the SSA22 proto-cluster region. The intrinsic relation 
should be determined for the general LAE population at 
lower z, where it is irrelevant to the reionization. 

The possible implications for reionization described 
above are still speculative because there are many un- 
known factors. First, we did not take into account the in- 
herent galactic evolution of LAEs. It is suggested by pre- 
vious studies that the number density of LAEs between 
z = 3 and 6 is almost unchanged, which is in contrast to 
the LBG population, w hich undergoes significant LF evo- 
lution from z = 3 to 6 llOuchi et alJl2004t iBouwens et alJ 



l2005UYoshida et al J 1200^1 . However, less than 20 LAEs 
have been confirmed spectroscopically, even at z ~ 3; 
this is not eno ugh to determine the LF accurately. On 
the other hand. IBouwens et al.l (2005) have found an ev- 
idence for evolution of the rest-frame continuum UV LF 
between z ~ 3 and 6. Also, there is no consensus on the 
spatial distribution of LAE at low redshifts. They might 
have an intrinsically homogeneous distribution. The in- 
herent evolutions of Lya flux, dust content, and neutral 
gas fraction inside LAEs themselves would also decrease 
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the number density from z — 5.7 to 6.5. The complex es- 
cape mechanism of ionizing radiation from galaxies is un- 
clear and depends strongly on assumed parameters, such 
as the escape fraction of Lya photons, amount of dust, 
galactic wind, and star formation activity. Our observed 
z6p5LAE LF could be more cons istent with the qu i escent 
hierarchical model prediction bv lLe Delliou et all l|2005t) 
than with their model including an ongoing starburst. 
Finally, the IGM physical conditions during the reioniza- 
tion epoch are unknown, and no observational evidence 
has been found for the cosmological H n region, which is 
the fundamental prediction of today's reionization mod- 
els. Although there are many ingredients to be consid- 
ered, we conclude that our conjectures about reionization 
are plausible and provide a reasonable explanation of our 
results. 

8. CONCLUSIONS 

We carried out spectroscopic observations with Subaru 
and Keck to identify z6p5LAEs that were selected by NB 
excess at ~ 920 nm. Our conclusions can be summarized 
as follows. 

1. We have identified eight new z6p5LAEs based 
on their significantly asymmetric Lya emission profiles. 
This increases the sample of spectroscopically confirmed 
z = 6.5 LAEs in the SDF to 17. 

2. We have constructed a large, homogeneous spectro- 
scopic sample from the photometric sample of 58 LAE 
candidates to determine the Lya LF at z = 6.5. The LF 
reveals an apparent deficit, at least at the bright end, 
compared to that at z — 5.7. The L* difference between 
z = 6.5 and 5.7 is ~ 0.75 mag for fixed a = —1.5. The 
difference has 3 a significance, which decreases to 2 a 
when we take into account cosmic variance. There may 
also be a decrease in comoving number density of LAEs 
from z = 5.7 to 6.5 at the faint end, although this conclu- 
sion could be modified by further follow-up spectroscopy. 

3. The rest-UV continuum LF of our LAE sample at 
z = 6.5 is almost the same as that of the LAE sample at 
z = 5.7 and the i-dropout sample at z ~ 6, even at their 
bright ends. 

4. The spatial distribution of z6p5LAEs was found to 
be homogeneous over the field, based on three indepen- 
dent methods to quantify the clustering strength. We 
cannot rule out the possibility that we happen to see a 
very homogeneous region of the z = 6.5 universe. 

5. The composite spectrum of the 12 z6p5LAE objects 
with high spectral resolution clearly reveals an asymmet- 
ric Lya profile with an extended red wing. The profile 
can be explained by either a galactic wind model com- 
posed of double Gaussian profiles or by a reionization 
model expected for a typical profile during the reioniza- 
tion epoch. 



6. Our results could imply that the reionization of the 
universe has not been completed at z = 6.5. The decline 
of the z6p5LAE LF implies < x^ M < .45 based on 
the IGM dynamical model of lSantoa l)2004[) . We conjec- 
ture that we are observing the final stage of reionization 
at z = 6.5, when the LAE overdense regions were still ob- 
scured by surrounding thick, neutral IGM clouds, which 
qualitatively agrees with our results of deficient bright 
LAEs and their homogeneous spatial distribution. 

Our z6p5LAE spectroscopic sample is not yet large 
enough to make a more precise comparison with LFs at 
lower redshifts. At the moment, it is not clear whether 
the true LF of z6p5LAE is closer to our upper or lower 
limit estimates. Thus, we do not conclude whether dif- 
ferences in L* or <fi* dominate in the disagreement of LF 
between z = 6.5 and 5.7. In this study, we could not 
constrain a, the faint-end slope of the LF. The faint end 
of the z6p5LAE sample would critically determine the 
LAE contribution as ionizing sources to the photon bud- 
get of cosmic reionization. Even if the number density 
of LAEs is low, a significant number of another star- 
forming population, LBGs, if they exist at this high-z 
epoch, could complete the reionization. Alternatively, 
only a small number of galaxies with huge star forma- 
tion rates at very h igh redshift can rcionizc the universe 
(Panada et al. 2005; Mobasher et al. 2005). It is too dif- 
ficult to sample LBGs at exactly the same redshift as 
LAEs, whereas it is not certain what fraction of LBGs 
at this epoch shows strong enough Lya emission to be 
observed as LAEs. Understanding the evolutionary con- 
necti on between LBGs and L AEs is linked to this prob- 
lem. Shimasaku et al. (2006) found from the rest UV 
luminosity LF of the z = 5.7 LAE sample that ~ 80% 
or more of the LBG population would have strong Lya 
emission at z ~ 6. This study shows that the rest-UV 
luminosity LF at z — 6.5 agrees with that of z — 5.7, sug- 
gesting the same high fraction. The faint end of the LF 
of z6p5LAEs may barely be determined with the spectro- 
scopic capability of today's 8 m telescopes. Nevertheless, 
steady efforts toward further spectroscopic confirmation 
are certainly required for z6p5LAEs. 
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APPENDIX 



SKEWNESS: THE ASYMMETRY INDICATOR OF HIGH-Z LYa EMISSION LINES 

We describe here a statistic skewness S to measure the asymmetry of high-z Lya emission lines. This model- 
independent indicator is free from fitting procedures that sometimes require smoothing of noisy spectra. 

The expression S is a popular statis tic, defined as the third moment of a distribution function, which describes its 
asymmetry (see also lKurk et al.ll2004() . Here, we regard the observed spectrum, which is basically a two-dimensional 
array of the flux (/,) and the pixel (xj), as a distribution function with an array size of n. Then S is defined as 
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respectively. ^^^^_^^^^ ^_^ 

This statistic is free from fitting procedures such as the a\ and aj presented bv lRhoads et alJ l)2003l) . Their asymme- 
try estimation is based on two-component Gaussian profile fitting for the red and blue sides of the emission. However, 
it is sometimes too difficult to determine accurately the peak wavelength X p of the line, or Aio, where the flux drops 
to 10% of its peak value, since they strongly depend on the resolving power and quality of the data. 

We now estimate the error in S, which can be analytically derived. We here assume Sxi ~ and approximately 
regard the lst-order derivative of skewness as its error, 
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We can assume that the flux error does not depend strongly on the wavelength, as in the case of our NB filter 
coverage, where the night-sky window is almost free from OH emission lines. Under this assumption, the 5 ft can be 
regarded as Sfi ~ Sf = const, dominated by readout noise around the emission line in the spectrum. 

This error estimate is confirmed to agree well with the rms fluctuation of S evaluated by a Monte Carlo realization, 
in which the line-profile model was disturbed with random artificial errors as large as Sf. Our estimate of skewness 
error in equation (A4) gives a useful analytic formula, although a more strict error estimate can only be achieved 
directly with Monte Carlo simulations. ^^^_^^^^ ^_^ 

We compared this statistic S with a\ : proposed by Rhoad set al.l l)2003f ). for our 7Vi?921-excess sample in Figure lAd2l 
We calculated S in the effective wavelength range of an emission line, where ft has an apparent positive signal above 
the continuum (sky) level. The slight change in the effective wavelength range does not significantly affect the result. 
Foreground emitters shown in Figure IA12I were definitely identified by their multiple spectral lines. That is, an Ha 
emitter has corresponding [O in] doublets around 7016 A, and [O in] and [O II] emitters show apparent doubles by 
themselves (and sometimes H(3 for the [O in] case). As expected, [O n], [O in], and Ha emitters are distributed around 
S = and a w = 1. The resolved [O n] doublet lines are expected to show negative S because A3726 is typically weaker 
than A3729. There is a population (filled circles) that have actually larger positive S and a w , indicating statistically 
asymmetric lines with a broad red wing. We recognize these as Lya emitters that have S > 0.15. Almost all the 
single-line emitters (triangles) show relatively low S, which indicates that they are likely to be [O n] emitters. The 
large scatter of S for these single-line emitters was caused by the low S/N of their spectra. 

To find a further adequate indicator sensitive to an asymmetry, we introduce weighted skewness S w , which combines 
the S indicator and Rhoads et al. (2003) 's method. High-z Lya emission usually has a wider FWHM in the observed 
frame than that of foreground emitters. We define the weighted skewness S w as 



S w — S'(Aio,r — Aio,b), 



(A7) 



where Aio, r and Aio,& are the wavelengths where the flux drops to 10% of its peak value at the redder and bluer 
sides of the Lya emission, respectively. The asymmetric index S is a dimensionless quantity, whereas the S w has a 
dimension of wavelength (here measured in angstroms) . Figure IA13I shows the correlation between S w and a w , where 
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Fig. A12. — Top: Comparison of the asymmetry indicator S with a\ proposed bv IRhoads et all (2003). The filled circles denote 
objects classified as z6p5LAEs in this study, whereas squares, triangles, and open circles, are the apparent [O III], [O o], and Ho emitters, 
respectively. The crosses denote single-line emitters. Bottom: The S distribution for z6p5LAEs (right histogram) and foreground emitters 
(left histogram). [See the electronic edition of the Journal for a color version of this figure.] 
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Fig. A13. — Same as Fig. IA12I but for S w - The errors in Aio, r and Aio,f, are not taken into account, and the error in S w is simply 
derived using <5S(Aio. r — Aiq (,). [See the electronic edition of the Journal for a color version of this figure.] 




the symbols are the same as those in Figure IA12I We found that all the apparent foreground emitters have S w < 3. 
The S w can distinguish between Lya and other lines more effectively than S. We can set the critical value S w = 3 to 
distinguish LAEs and foreground emitters, although there could be more or less contamination of LAEs at S w < 3. 
This is a conservative critical value for LAEs in the sense that it assumes low contamination and low completeness. 

In summary, our proposed statistic S is free from fitting procedures and has an analytical error estimation formula. 
In addition, we made use of a revised indicator S w that proves to be more powerful when combined with the line width 



Lya emitters at z = 6.5 
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determined from Gaussian fitting. 
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